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Features  of  light  absorption  are  critical  to  optical  aspects  of  water  quality  and  in  regulating  the  signal  available  for 
remote  sensing.  Spectral  characteristics  and  spatial  patterns  of  light-absorbing  components,  and  their 
relationships  with  optically  active  constituents,  are  documented  for  the  Sturgeon  River,  Keweenaw  Bay.  and 
Lake  Superior  based  on  analyses  of  samples  collected  on  two  cruises  (2006  and  2007, 20  sites).  The  absorption 
coefficient,  a  (m“  ^ ),  is  partitioned  according  to  the  additive  components  (Ox)  of  colored  dissolved  organic  matter 
(ocdom).  non-algal  particles  (aNAp).  phytoplankton  ((4).  and  water  itself  (awi  known).  The  role  of  minerogenic 
particles  and  their  iron  content  in  regulating  Onap  is  evaluated  based  on  paired  measurements  by  an  individual 
particle  analysis  technique  (Peng  et  al..  2009),  through  empirical  analyses  and  Mie  theory  calculations  of 
absorption  by  these  particles  (Om).  Spectral  characteristics  of  Qnap  and  were  consistent  with  those  reported  for 
other  case  2  (i.e.,  phytoplankton  not  dominant)  systems.  However,  the  slope  values  that  describe  Ocdom  spectra 
for  the  bay  and  the  lake  were  unusually  low,  suggesting  an  atypical  composition  for  the  lake's  CDOM.  The 
dominant  absorbing  component  in  the  blue  wavelengths  was  CDOM.  representing  ^75%  of  a  at  a  wavelength  of 
440  nm  at  all  sites  in  the  2006  survey.  A  general  gradient  in  both  Ocdom  and  Onap  extended  from  the  Sturgeon 
River,  through  the  bay,  into  eastern  Lake  Superior  in  that  survey.  Relationships  between  and  optically  active 
constituents  were  within  the  broad  ranges  reported  for  other  case  2  systems.  Minerogenic  particles,  related  to 
their  iron  content,  are  demonstrated  to  be  an  important  component  of  Onap- 
©  2010  International  Association  for  Great  Lakes  Research.  Published  by  Elsevier  B.V.  All  rights  reserved. 


Introduction 

The  magnitude  of  the  light  absorption  coefficient  (a;  unit  of  1/ 
length.  m~^),  an  inherent  optical  property  (lOP;  i.e.,  independent  of 
the  geometry  of  the  light  field),  is  the  fraction  of  the  radiant  flux 
absorbed  divided  by  the  thickness  of  the  layer  (Kirk,  1994).  The 
magnitude  of  a  varies  as  a  function  of  wavelength  (\),  according  to  the 
composition  of  water  with  respect  to  light  absorbing  constituents 
(Kirk,  1994;  Babin  et  al..  2003b).  The  magnitude  and  spectral 
character  of  a  are  determinants  of  apparent  optical  properties 
(AOPs;  those  that  depend  on  the  geometry  of  the  light  field),  such 
as  the  diffuse  attenuation  coefficient  and  remote  sensing  reflectance. 

Usually,  a(\)  is  partitioned  into  four  additive  components  (Ox; 
Babin  et  al..  2003b) 

a(\)  =  a^(\)  +  a4,(\)  -f  aNAp(^)  +  Qcdom(^)  (1) 

that  include  contributions  from  pure  water  (w),  phytoplankton  (<J)), 
non-algal  particles  (NAP),  and  colored  dissolved  organic  matter  (CDOM. 
or  gelbstoff).  Absorption  by  pure  water  increases  in  red  wavelengths, 
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and  is  well  defined  (Pope  and  Fry,  1997).  Phytoplankton  demonstrate  a 
broad  absorption  peak  at  -440  nm  and  a  more  well  defined  maximum 
at  676  nm  (Kirk,  1994;  Babin  et  al.,  2003b).  Both  Onap  and  Ocdom 
decrease  exponentially  with  increasing  \  over  the  visible  wavelengths 
(Babin  et  al..  2003b;  Perkins  et  al.,  2010).  Knowledge  of  a,  its 
components  (Ox).  and  their  spectral  characteristics  is  important  for 
developing  and  validating  remote  sensing  algorithms  for  water  quality 
and  mechanistic  optics  models  for  the  underwater  light  field  (Effler  et 
al..  2008),  and  advancing  related  initiatives  for  the  Laurentian  Great 
Lakes  (Cons  and  Auer,  2004;  Cons  et  al..  2008).  Establishment  of 
relationships  between  Ox  and  concentrations  of  optically  active 
constituents,  such  as  total  suspended  solids  (TSS)  and  chlorophyll  a 
(Chi),  is  necessary  to  support  mechanistic  approaches  to  estimate  their 
patterns  through  remote  sensing. 

Case  1  waters  (open  oceans)  are  those  where  patterns  in  optical 
properties  are  driven  by  distributions  of  phytoplankton,  and  NAP  and 
CDOM  covary  with,  and  are  subordinate  to.  the  phytoplankton 
component  (Morel  and  Prieur,  1977).  In  contrast,  in  case  2  waters, 
such  as  coastal  marine  and  lacustrine  systems.  NAP  and  CDOM  do  not 
necessarily  covary  with  phytoplankton  and  may  not  be  subordinate. 
Direct  quantification  of  Ox  in  lacustrine  systems  remains  rare  (Binding 
et  al.,  2008;  Perkins  et  al..  2009,  2010),  despite  the  establishment  of 
standardized  laboratory  protocols  for  a^.,  Onap  ^nd  Ucdom  for  marine 


0380-1330/$  -  see  front  matter  C  2010  International  Association  for  Great  Lakes  Research.  Published  by  Elsevier  B.V.  All  rights  reserved, 
doi:  1 0. 1 0 1 6/j.jglr.20 1 0.08.00 1 


S.W.  Effler  et  al.  /  Journal  of  Great  Lakes  Research  36  (2010)  656-665 


657 


initiatives  (Mitchell  et  al.,  2003).  Resolution  and  quantification  of  NAP 
have  received  the  least  attention  of  the  components  (Babin  et  al., 
2003b).  Moreover,  there  exists  great  uncertainty  concerning  NAP's 
composition,  particularly  related  to  contributions  of  minerogenic 
versus  organic  (e.g.,  detritus)  particles  (Bowers  et  al.,  1996;  Babin  et 
al.,  2003b;  Binding  et  al.,  2008).  Babin  and  Stramski  (2004)  reported 
that  absorption  by  minerogenic  particle  populations  (a^)  was  largely 
attributable  to  the  content  of  iron  (Fe).  Recently,  advancements  have 
been  made  in  estimating  the  minerogenic  scattering  coefficient  (bm) 
based  on  an  individual  particle  analysis  technique  that  provides  both 
morphological  and  chemical  characterizations  (Scanning  electron 
microscopy  interfaced  with  Automated  image  and  X-ray  analyses 
(SAX)),  coupled  with  forward  Mie  theory-based  calculations  (Peng 
and  Effler.  2007;  Peng  et  al..  2007,  2009).  The  integrated  SAX-Mie 
theory  approach  also  represents  a  basis  to  directly  estimate  Um, 
mediated  through  the  imaginary  (absorption)  part  of  the  refractive 
index  (Wozniak  and  Stramski,  2004). 

Lake  Superior  has  been  described  as  an  optically  complex  system 
that  represents  a  challenge  for  robust  remote  sensing  of  water  quality 
(Li  et  al..  2004).  Partitioning  of  a  into  Ux  contributions  for  visible 
wavelengths  for  multiple  sites  has  not  been  reported  for  this  lake. 
Minerogenic  particles,  mostly  clay  minerals  in  the  1  -1 0  pm  size  range, 
were  found  to  be  an  important  component  of  backscattering  in 
eastern  Lake  Superior  and  Keweenaw  Bay  in  July  of  2006  (Peng  et  al.. 
2009).  Rigorous  partitioning  of  a  according  to  Ox  and  establishment  of 
relationships  between  and  concentrations  of  optically  active 
constituents  have  been  recommended  for  the  lake  to  support  remote 
sensing  initiatives  (Cons  and  Auer,  2004;  Cons  et  al..  2008).  This  paper 
documents  spectral  characteristics,  spatial  patterns,  and  relationships 
with  optically  active  constituents  for  for  Lake  Superior  and 
Keweenaw  Bay  based  on  surveys  in  2006  and  2007.  The  results  are 
considered  in  the  context  of  the  modest  literature  that  describes  for 
lacustrine  systems  and  other  case  2  waters,  and  remote  sensing  issues 
for  the  lake.  Symbols  and  abbreviations  adopted  in  the  paper  (Table  1 ) 
are  consistent  with  those  used  by  the  optics  research  community. 
Particular  attention  is  given  to  the  role  of  minerogenic  particles  and 
their  Fe  content  in  influencing  Qnap.  based  on  paired  SAX  results  (Peng 
et  al.,  2009)  and  Mie  theory  calculations. 

Methods 

Study  sites 

This  study  is  based  on  analysis  of  samples  collected  in  two  surveys 
conducted  in  July  (5-14)  2006  and  August  (18-24)  2007.  The  2006 
survey  included  1 1  sites  extending  from  the  mouth  of  the  Sturgeon 
River  in  the  Portage  Lake  portion  of  the  Keweenaw  Waterway  (KWO), 
down  to  the  waterway's  exit  into  Keweenaw  Bay  (KWl),  along  the 
primary  axis  of  the  bay  (KBl  -KB6),  and  out  into  eastern  Lake  Superior 
(three  sites,  ELS1-ELS3;  Fig.  1).  Site  designations  (except  KWO)  are 
those  used  by  Peng  et  al.  (2009).  The  bay  is  36  km  long,  with  widths 
ranging  from  3  to  1 8  km.  The  bay  sites  represent  a  subset  of  the  twelve 
included  in  a  study  of  its  intrinsic  water  color  (Cons  and  Auer,  2004). 
The  Sturgeon  River  has  been  described  as  humic  (Ma  and  Green, 
2004)  and  highly  stained  (Cons  and  Auer,  2004).  This  river  is  the 
largest  tributary  entering  the  bay,  with  an  average  annual  flow  rate  of 
26  m^s'^ 

The  2007  survey  included  nine  sites  (Fig.  1 )  that  represent  a  subset 
of  those  monitored  in  USEPA's  long-term  program  (GLNPO,  2009). 
Site  designations  for  the  2007  samples  correspond  to  those  used  in 
that  monitoring  program. 

Measurements 

Near-surface  (within  upper  2  m)  samples  were  collected  and 
processed  the  day  of  collection  for  subsequent  laboratory  analyses. 


Table  1 

Abbreviations  and  symbols. 


Abbreviations 

AOPs 

Apparent  optical  properties:  those  that  depend  on  the 
underwater  light  held 

CDOM 

Colored  dissolved  organic  material 

Chi 

Chlorophyll  a 

DOC 

Dissolved  organic  carbon 

FeXRI 

Iron  X-ray  relative  intensity  (%) 

FSS 

Fixed,  or  inorganic,  suspended  solids 

NAP 

Non-algal  particles 

PAV^ 

Projected  area  concentration  of  minerogenic  particles  (m~  *) 

SAX 

Scanning  electron  micro.scopy  interfaced  with  automated 
image  and  X-ray  analyses 

TSS 

Total  suspended  solids 

VSS 

Volatile  (organic)  suspended  solids 

XRI 

X-ray  relative  intensity  (%) 

Symbols 

a(\) 

Total  spectral  absorption  coefficient  (m“  '),  wavelength  (A) 
sometimes  omitted 

ax(X) 

Spectral  absorption  coefficient  for  component  “x"  (m“  ’);  “x" 

(as  subscript)  can  be  CDOM.  NAP.  detritus  (d).  particles  (p). 
minerogenic  particles  (m).  phytoplankton  (tp).  and  water  (w) 

flxCXf) 

a  for  component  “x"  (CDOM  or  NAP)  at  reference 
wavelength  A^  (m“ ') 

flax)M.Doc(440) 

DOC-specific  Ocdom  at  A  =  440  nm  (m^  g*  ’) 

aNAP.Tss(440) 

TSS-specific  Onap  at  A  =  440  nm  (m^g“’) 

aNAP.Fss(440) 

FSS-specific  Anap  at  A  =  440  nm  (m^  g“  ’ ) 

aJ.T3s(440) 

TSS-specific  flp  at  A  =  440  nm  (m^  g“  ’) 

<.chi(440) 

Chl-specific  a,p  at  A  =  440  nm  (m^g“  ') 

<chi(676) 

Chl-specific  a^  at  A  =  676  nm  (m^  g“ ') 

bm 

Scattering  coefficient  for  minerogenic  particles  (m~  ’) 

dj 

Size  of  minerogenic  particle  j  (pm) 

m 

Relative  (to  water)  complex  refractive  index 

n 

Real  part  of  m 

n' 

Imaginary  part  of  m 

Nm 

Number  of  minerogenic  particles  per  unit  volume  of  water 
(numberm~^) 

PAmo 

Projected  area  of  minerogenic  particle  j 

Qamj 

Absorption  efficiency  factor  for  minerogenic  particle  J 

Scdom 

Slope  of  the  acDOM(A)  spectrum  (nm“  ’) 

SrjAP 

Slope  of  the  aNAp(A)  spectrum  (nm  “  ’) 

V 

Water  sample  volume  for  SAX  (m“  ^) 

K 

Reference  wavelength  of  light  (nm) 

Measurements  of  concentrations  of  optically  active  substances  were 
conducted  on  sample  aliquots  for  total  (TSS;  Whatman  934-AH, 
1.5  pm  pore  size)  and  fixed  (FSS,  after  550  “C;  i.e.,  inorganic) 
suspended  solids  (Clesceri  et  al..  1998),  dissolved  organic  carbon 
(DOC)  (Clesceri  et  al.,  1998;  for  2007  samples  only),  and  chlorophyll  a 
(Chi;  Parsons  et  al..  1984). 

The  methods  for  optical  measurements  generally  tracked  those 
specified  for  satellite  ocean  color  calibration/validation  activities 
(Mitchell  et  al.,  2003).  Determinations  of  Ocdom.  Qnap.  and  04,  were 
made  on  samples  from  2006,  while  measurements  were  limited  to 
flcDOM  in  2007.  Lake  water  was  filtered  under  low  vacuum  on  pre¬ 
rinsed  (deionized  water)  0.2  pm  pore  size  Millipore  membranes  for 
CDOM  absorption  measurements  and  refrigerated  (4  "C)  in  the  dark 
until  analyzed.  Absorbance  of  the  filtered  water  was  measured  in  a 
10  cm  quartz  cuvette  from  400  to  700  nm  at  2  nm  increments  with  a 
dual  beam  spectrophotometer  (Perkin  Elmer,  Lambda  40),  with 
deionized  water  used  as  the  reference.  A  baseline  correction  was 
made  by  subtracting  the  absorbance  value  averaged  over  a  14nm 
interval  around  693  nm  from  all  the  spectral  values.  These  spectra 
were  described  according  to 

^*CDOm(^)  =  OcDOm(\)^^  ^cdom(>^  ^))  ^2) 

where  Ar  (nm)  is  the  reference  wavelength,  Ocdom  (K)  is  the  CDOM 
absorption  at  Ap.  and  Scdom  (nm“M  is  the  slope  of  the  Ucdom  (X) 
spectrum.  Values  of  Scdom  were  estimated  through  a  nonlinear 
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Fig.  1.  Map  of  Lake  Superior  with  sampling  sites  in  the  Sturgeon  River  (KWO)  and  its 
inflow  zone  (KWl),  Keweenaw  Bay  (KB.  6  sites),  eastern  Lake  Superior  (ELS.  3  sites. 
2006).  and  lake-wide  (SU.  9  sites.  2007).  Site  monitored  by  Li  et  al.  (2004)  is  also 
shown. 

regression  fitting  approach  (Twardowski  et  al.,  2004)  for  the 
exponential  function  (Eq.  (2),  with  \r  =  440nm)  over  the  400  to 
500  nm  wavelength  interval  (Perkins  et  al..  2(X)9.  2010).  It  should  be 
noted  that  this  fitting  method  differs  from  that  of  taking  the  logarithm 
of  flcDOM  absorption  and  obtaining  Scdom  ^s  the  slope  from  a  linear 
least-squares  fit  (Twardowski  et  al.,  2004). 

Water  samples  were  filtered  onto  25  mm  diameter  glass  fiber 
filters  (Whatman  GF/F,  nominal  pore  size  of  0.7  pm)  at  low  vacuum 
for  particle  absorption  measurements.  Filters  were  stored  in  opaque 
containers  in  the  dark  and  refrigerated  until  analyzed.  The  absorption 
spectrum  of  particles  retained  on  the  filters.  ap(A),  was  determined 
using  the  method  of  Lohrenz  (2000).  Particle  absorbance  was 
measured  between  400  and  750  nm  at  2  nm  increments  using  a 
dual  beam  spectrophotometer  equipped  with  a  100  mm  integrating 
sphere  (Perkins  Elmer,  Lambda  18).  Correction  for  path  length 
amplification  was  made  using  the  relationship  developed  by  Cleve¬ 
land  and  Weidemann  (1993).  aNAp(^)  was  measured  after  pigment 
bleaching  with  sodium  hypochlorite  (Ferrari  and  Tassan,  1999).  The 
final  determinations  of  ap(\)  and  UnapC^)  were  made  by  subtracting 
the  average  of  measured  values  over  the  740  to  750  nm  range  from 
the  spectral  observations.  The  phytoplankton  absorption  spectrum, 
(\),  was  obtained  through  a  residual  calculation 

=  ap(K)-aNAp(M  (3) 

This  determination  of  a^{\)  includes  absorption  by  all  pigments 
bleached  by  sodium  hypochlorite.  The  Onap  spectra  were  described  in 
an  analogous  manner  to  Ucdom  according  to 

(4) 


where  Unap  (V)  is  NAP  absorption  at  \r  (440  nm),  and  Snap  is  the  slope 
of  the  Onap  (^)  spectrum.  Values  of  Snap  were  estimated  through  the 
same  nonlinear  regression  fitting  approach  used  for  Scdom  for  the 
exponential  function  (Eq.  (4),  with  \r  =  440  nm)  for  data  from  482  to 
618  nm  and  712  to  750  nm  (Babin  et  al.,  2003b).  The  exclusion  of 
400-480  and  620-710  nm  ranges  avoids  the  potential  effects  of 
residual  pigment  absorption  that  may  have  persisted  after  the  sodium 
hypochlorite  treatment  (Babin  et  al..  2003b). 

SAX  data  and  absorption  by  minerogenic  particles 

Results  from  SAX  characterizations  of  the  2006  samples,  in  the 
context  of  light-scattering  implications,  were  reported  by  Peng  et  al. 
(2009).  The  details  of  sample  preparation  and  instrumentation 
methodology  for  SAX  have  been  described  elsewhere  (Peng  and 
Effler.  2005,  2007).  Approximately  2000  individual  particles  were 
analyzed  for  each  sample  for  both  composition  and  morphology 
characteristics  (Peng  et  al..  2009).  SAX  provided  information  on 
chemical  composition,  number  concentration,  and  size  distribution  of 
minerogenic  particles  for  forward  estimates  of  both  bm  ^nd  0^. 

Composition  of  individual  minerogenic  particles  is  represented  by 
the  X-ray  relative  intensities  (XRl,  in  percentage)  of  16  elements 
[including  aluminum  (Al),  silica  (Si),  calcium  (Ca).  and  iron  (Fe)|. 
Based  on  composition,  particles  were  classified  into  five  predefined 
generic  types:  clay,  quartz,  silica  (Si)-rich.  calcium  (Ca)-rich.  and 
miscellaneous  (Peng  et  al.,  2009).  The  projected  area  of  minerogenic 
particles  per  unit  volume  of  sample  (PAV^),  determined  with  SAX 
(Peng  and  Effler.  2007),  is  a  fundamental  attribute  of  the  light¬ 
attenuating  properties  of  these  particles  (Peng  and  Effler.  2007;  Peng 
et  al.,  2009).  The  minerogenic  particle  populations  of  the  2006  survey 
were  dominated  by  clay  minerals,  representing  on  average  72%  of  the 
total  PAVm  at  each  of  the  sites  (Table  2).  Such  particles  have 
terrigenous  origins  (Peng  et  al.,  2009). 

The  SAX-Mie  theory  calculation  approach  was  used  here  to  make 
forward  estimates  of  based  on  the  same  individual  particle 
characterizations,  according  to 

an.(M  =  (5) 

where  V  is  the  water  sample  volume.  is  the  number  of  minerogenic 
particles  per  unit  volume  of  water,  (iimj  is  the  absorption  efficiency 
f^actor  (determined  from  Mie  theory  calculations;  Bohren  and  Huff¬ 
man.  1983)  for  minerogenic  particle  j,  n'  is  the  imaginary  portion  of 
the  complex  refractive  index  (m  =  n  — in',  where  n  the  is  the  real 
portion),  \  is  the  wavelength,  dj  is  the  size  of  particle  j,  and  PAmj  is  the 
projected  area  of  minerogenic  particle  j.  Values  of  n'  are  more 
uncertain  than  those  for  n  (Wozniak  and  Stramski,  2004).  The 
relationship  between  minerogenic  particles  and  Onap  was  investigated 
in  two  ways:  (1)  empirically,  through  regression  analyses  between 
□nap  3nd  PAVm  ^nd  Qnap  ^nd  average  FeXRl  of  minerogenic  particles 
(not  reported  by  Peng  et  al.  (2009)),  and  (2)  mechanistically,  through 
a  series  of  Mie  theory  calculations  of  0^(440)  based  on  the  paired  SAX 
data  (Peng  et  al..  2009).  Various  cases  of  n'  were  adopted  in  the  Mie 
theory  calculations  that  included  the  approach  to  equivalence  (i.e.. 


Table  2 

Summary  of  average  minerogenic  particle  chemistry  according  to  five  classes,  in  the 
context  of  contributions  to  PAV^,. 


Chemical  class 

Average  contribution  (%) 

n' 

Clay 

72.3 

4  scenarios'* 

Quartz 

7.7 

0.0002 

Si-rich 

10.2 

4  scenarios 

Ca-rich 

3.0 

0.0002 

Miscellaneous 

6.9 

4  scenarios 

*  n'  =  0.003,  0.006,  0.008,  0.010. 


aNAp(X)  =aNAp(X,)e'-^'’'-’^'’ 
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closure)  with  Unap-  The  adopted  n'  values  were  guided  by  the  related, 
albeit  limited,  literature  (Patterson  et  al.,  1977;  Wozniak  and 
Stramski,  2004).  Values  of  n'  were  assumed  to  be  equal  for  the  clay, 
Si-rich,  and  miscellaneous  particle  classes  for  the  various  cases  of  n\ 
while  the  quartz  and  Ca-rich  classes  had  much  lower  values  that  were 
held  constant  (Table  2;  Wozniak  and  Stramski.  2004). 

Results 

Spectral  characteristics 

Representative  spectra  for  the  three  measured  components  (a^) 
are  illustrated  for  site  KB3  (Fig.  2,  note  scale  changes  for  Ux). 
Exponential  decreases  with  increasing  wavelength  were  observed 
for  OcDOM  (Pig-  2a)  and  Onap  (Pig-  2b).  The  non-linear  regression  fits 
generally  performed  well  in  representing  these  spectra  over  the 
designated  calculation  wavelength  intervals  for  the  slope  values  (Fig. 
2a  and  b).  All  phytoplankton  spectra  demonstrated  two  maxima,  a 
broad  one  at  -440  nm  and  the  other  better  defined  at  676  nm 
(Fig.  2c).  Instrument  noise  was  apparent  for  both  Qnap  (Pig-  2b)  and  a^ 
(Fig.  2c)  spectra  because  of  the  relatively  low  absorption  signals  for 
these  components  associated  with  their  low  concentrations. 

The  Scdom  values  for  KWO  and  KWl  were  the  highest  of  the  2006 
survey,  though  these  were  approached  by  the  observation  for  KBl 
(Table  3).  Values  for  the  lake  and  bay  sites  were  in  the  range  of 
0.01 1 1-0.0160  nm~  ^  the  mean  for  these  sites  was  0.0126  nm~  \  the 
coefficient  of  variation  (cv)  was  1 1  %.  The  mean  Snap  for  these  sites  was 
similar.  0.0131  nm'^  (cv  =  13%).  Significant  relationships  between 
5cdom  dnd  OcDOM.  ^nd  Snap  and  Unap.  were  not  observed  for  the  lake 
and  bay  sites. 

Values  of  Scdom  for  the  lake  sites  of  the  2007  survey  (Table  4)  were 
generally  similar  to  those  of  the  KB  and  ELS  sites  of  2006  (Table  3).  The 
2007  observations  were  in  the  range  of  0.0089  (SU23)  to  0.0139  nm“  ^ 
(SU22)  for  all  sites,  with  an  average  of  0.011  nm”^  (cv=14%).  a^ 


Fig.  2.  Representative  spectra  for  light  absorbing  components  for  site  KB3:  (a)  Qcdom. 
with  fit,  (b)  Qnap.  with  fit,  and  (c)  a^,.  Note  scaling  difference  for  a*. 


spectra  were  characterized  by  the  ratio  of  the  blue  (\  =  440  nm)  to  red 
(\  =  676nm)  peaks  (Table  3,  e.g.,  Babin  et  al..  2003b:  Perkins  et  al.. 
201 0).  All  but  one  of  the  ratio  values  were  in  the  range  of  2.06-3.1 9,  with 
an  average  of  2.54  (cv  =  23%)  excluding  site  KB3  (4.24). 

Absorption  by  components 

Spatial  patterns  of  Ux  are  described  for  A  =  440  nm,  the  commonly 
adopted  blue  reference  wavelength  (Kirk.  1994).  Values  of  Ocdom 
(440)  decreased  in  a  generally  progressive  manner  from  the  Sturgeon 
River  to  the  open  waters  of  the  lake  in  the  2006  survey,  according  to 
river > inflow  mixing  zone > bay > eastern  Lake  Superior  (Fig.  3a,  note 
scale  change  for  aj,  though  substantial  spatial  variability  was 
manifested  (Table  5).  On  average,  the  KB  site  levels  were  70%  higher 
than  the  ELS  sites.  Though  aNAp(440)  levels  were  substantially  lower,  a 
similar  spatial  pattern  was  observed  (Table  5,  Fig.  3b).  The  river  Onap 
(440)  value  exceeded  the  inflow  mixing  zone  level  by  more  than  2- 
fold,  which  in  turn  was  about  6-fold  greater  than  the  bay,  and  the  bay 
was  about  2-fold  higher  than  for  the  ELS  sites.  The  spatial  pattern  for 
a^  differed  in  that  average  values  for  the  KB  and  ELS  sites  were  nearly 
equal  (Fig.  3c).  The  spatial  pattern  of  a(440)  (Fig.  3d),  determined 
from  the  summation  of  Ux  (Eq.  (1 )),  generally  tracked  that  of  Ucdom 
(440)  (Fig.  3a),  the  dominant  absorbing  component. 

acDOM(440)  represented  >75%  of  a(440)  at  all  sites.  acDOM(440) 
was  80%  of  a(440)  for  KWO,  77%  for  KWl,  81%  for  the  average  of  the 
KB  sites,  and  75%  for  the  average  of  the  ELS  sites  ( Fig.  3e).  NAP  was  the 
second  largest  component  for  KWO  (12%)  and  KWl  (13%),  but  was 
less  than  a^  at  the  KB  (7  vs.  9%)  and  ELS  (9  vs.  13%)  sites.  aw(440)  was 
the  smallest  component  at  all  sites;  its  greatest  contribution  was 
where  a(440)  was  the  lowest  (3%,  on  average,  based  on  observations 
for  ELS  sites).  The  general  spatial  pattern  of  the  spectral  contributions 
of  Ox  is  illustrated  based  on  observations  for  KWO  (Fig.  4a),  KWl 
(Fig.  4b),  KB4  (Fig.  4c),  and  ELS2  (Fig.  4d).  The  fractional  contribution 
(Fig.  4)  format  serves  to  contrast  the  spatial  differences  as  well  as 
depict  challenges  faced  in  resolution  of  the  signals  from  the  various 
components,  particularly  a^.  a^  was  dominant  for  A>600nm.  This 
dominance  was  manifested  at  lower  wavelengths  as  the  other 
components  of  a  decreased.  The  contributions  Unap  ^nd  a^  were 
spectrally  relatively  featureless  at  KWO  and  KWl  (Fig.  4a  and  b) 
because  of  the  high  OcDOM-though  the  a^,  component  was  apparent  for 
KB4  and  ELS2  at  A<500  nm,  representing  about  12  and  20%  of  a  at 
these  maxima,  a^  was  larger  than  Qnap  ^nd  Ocdom  its  secondary 
peak  (A  =  676nm)  at  ail  sites,  but  its  contribution  to  a  at  this 
wavelength  was  small  compared  to  Ow.  The  contribution  of  a^  at 
A  =  676  nm  was  more  than  6-fold  greater  than  a<i,  at  KWO,  17-fold  at 
KWl.  53-fold  at  KB4.  and  40-fold  at  ELS2. 

Optically  active  constituents  and  relationships  with  a^ 

The  highest  Chi  in  the  2006  survey  was  in  the  Sturgeon  River 
(KWO)  (Table  3).  The  concentration  at  KWl  was  about  58%  lower,  but 
more  than  3-fold  higher  than  the  bay  concentrations.  Chi  levels  were 
higher  at  the  ELS  sites  than  those  at  the  KB  sites,  about  32%  on  average. 
TSS  and  FSS  concentrations  were  distinctly  higher  for  KWO  and  KWl 
than  the  bay  sites  (Table  3).  The  suspended  solids  data  for  the  bay 
sites  are  uncertain  because  of  the  low  weights  accumulated  on  the 
filters  in  the  analyses,  and  are  considered  applicable  only  on  a  bay- 
average  basis.  No  suspended  solids  measurements  were  available  for 
the  ELS  sites.  The  PAV^  data  are  those  reported  by  Peng  et  al.  (2009), 
except  for  KWO  (not  previously  reported).  A  distinct  gradient  in  this 
metric  of  minerogenic  particles  was  manifested,  according  to 
Sturgeon  River >  inflow  mixing  zone >  bay >  ELS  sites.  The  FeXRI  data 
presented  here  (Table  3)  from  SAX  were  not  reported  in  previous 
work.  The  highest  FeXRl  values  were  in  the  Sturgeon  River  and  its 
inflow  mixing  zone.  A  gradient  prevailed  in  the  bay,  with  decreases 
manifested  with  the  approach  to  the  open  waters.  The  levels  of  FeXRI 
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Table  3 

Optically  active  constituent  concentration  and  spectral  characteristics.  2006  survey. 


Chi  (mgm"’) 

7SS(gm-=) 

fSS  (gm-’) 

FeXRI 

FeXRI 

PAV„ 

ScDOM 

(nm"*) 

Snap 

(nm"*) 

a,>(440) 

04,(676) 

KWO 

3.85 

1.93 

0.93 

0.24 

8.8 

36 

0.0169 

0.0120 

2.083 

KWI 

1.63 

1.43 

1.03 

0.182 

9.4 

52 

0.0168 

0.0118 

2.920 

KB1 

0.48 

0.30 

0.10 

0.067 

6.3 

94 

0.0160 

0.0144 

2.540 

KB2 

0.54 

024 

0.04 

- 

- 

- 

0.0125 

0.0125 

3.187 

KB3 

0.42 

0.28 

0.07 

0.080 

6.2 

78 

0.0119 

0.0125 

4.243 

KB4 

0.45 

0.63 

0.33 

0.030 

5.3 

177 

0.0117 

0.0105 

2.852 

KBS 

0.48 

0.30 

0.12 

0.043 

4.1 

95 

0.01 1 1 

0.0131 

2.457 

KB6 

0.43 

0.37 

0.17 

0.050 

4.0 

80 

0.0126 

0.0165 

2.180 

ELS1 

0.69 

0.016 

4.1 

256 

0.0117 

0.0128 

2.550 

ELS2 

0.57 

0.026 

4.4 

169 

0.0127 

0.0121 

2.531 

ELS3 

0.81 

0.022 

4.0 

181 

0.0135 

0.0135 

2.063 

Mean 

0.94 

0.69 

0.35 

0.076 

5.7 

122 

0.0134 

0.0123 

2.691 

Std.  dev. 

1.02 

0.64 

0.40 

0.075 

2.0 

70.1 

0.002 

0.0016 

0.623 

were  nearly  equal  for  KBS,  KB6.  and  the  ELS  sites  (Table  3).  The 
relative  iron  content  of  the  minerogenic  particles  (e.g.,  on  a  projeaed 
area  basis,  indicated  by  the  ratio  FeXRl-^PAVmi  Table  3),  generally 
increased  from  the  river  to  the  bay,  and  from  the  bay  to  ELS  sites. 
Concentrations  of  DOC  were  in  the  narrow  range  of  1.3- 1.5  g  m“^  for 
the  2007  survey  of  Lake  Superior  (Table  4). 

Relationships  between  and  optically  active  constituents  were 
evaluated  as  “specific”  absorption  coefficients,  which  have  been 
normalized  by  the  constituent  concentration  of  interest  (e.g.,  Binding 
et  al..  2008).  These  relationships  support  comparison  to  values  of 
specific  absorption  coefficients  reported  for  other  systems.  The 
average  Chl-specific  a^p  for  \  =  440  (a*,}>.chi(440))  and  676  nm  (a%. 
chi(676))  for  the  2006  cruise,  excluding  KWO  and  KWl,  were  0.046 
and  0.017  m^mg“\  respectively:  the  corresponding  values  of  cv 
were  26  and  36%.  The  average  a*4,.chi(440)  for  the  KB  sites 
(0.050  m^mg~^)  was  significantly  (p<0.001)  higher  than  the  ELS 
sites  (0.037  m^  mg“M-  Average  specific  absorption  coefficients  for 
NAP  are  presented  based  on  both  TSS  and  FSS,  because  of:  ( 1 )  the 
ambiguity  associated  with  the  inclusion  of  phytoplankton  biomass 
within  the  TSS  measurement,  and  (2)  the  uncertainty  concerning  the 
relative  roles  of  minerogenic  versus  detrital  particles  in  regulating 
fliMAP.  The  average  specific  absorption  coefficient  values  for  NAP  based 
on  TSS  (a*NAP.TSs  =  aNAp(440)-^TSS)  for  KB  and  KW  sites  were  0.057 
(cv  =  57%)  and  0.10  m^  g~  \  respectively.  The  corresponding  specific 
absorption  coefficient  values  for  NAP  based  on  FSS  (a*NAP.  Fss  =  aNAP 
(440)  4-  FSS)  were  0.083  (cv  =  39%)  and  0.27  m^  g~  \  The  relationship 
between  Op  (a^y-hOnAp)  and  TSS  does  not  suffer  from  the  above 
ambiguities.  The  average  specific  absorption  coefficients  for  overall 
particle  absorption  (aV.Tss  =  ap(440)-^TSS)  for  the  KB  and  KW  sites 
were  0.13  (cv  =  36%)  and  0.17m^g“\  respectively.  The  average 
specific  absorption  coefficient  for  DOC  for  \  =  440nm  (qVdom.doc 


Table  4 

DOC  concentrations.  Scoom.  Qcdom  and  acDOM(440)  h-  DOC.  2007  survey. 


DOC 

ScDOM 

0cdom{AA0) 

acDOMiAAO) 

DOC 

(g  ni"3) 

(m-*) 

(m-*) 

(m2  g-') 

SUl 

1.3 

0.0108 

0.119 

0.092 

SU6 

1.5 

0.0107 

0.109 

0.073 

SU9 

1.3 

0.0103 

0.145 

0.112 

sun 

1.3 

0.0115 

0.132 

0.101 

SU12 

1.3 

0.0105 

0.150 

0.115 

SUM 

1.4 

0.0135 

0.073 

0.052 

SU20 

1.4 

0.0098 

0.208 

0.148 

SU22 

1.5 

0.0139 

0.157 

0.105 

SU23 

1.4 

0.0090 

0.234 

0.167 

Mean 

1.4 

0.0107 

0.145 

0.105 

Std.  dev. 

0.083 

0.0016 

0.049 

0.035 

(440)  =  acDOM(440)^DOC)  based  on  paired  measurements  for  the 
2007  cruise  was  0.10  m^g'^  (cv  =  37%). 

Minerogenic  particles,  Qnap,  and 

PAVm  was  a  strong  predictor  of  Qnap  for  the  KB  and  ELS  sites  of  the 
2006  survey,  explaining  69%  of  the  variability  according  to  linear  least 


Sites 

Fig,  3.  Spatial  patterns  of  ax(440)  for  2006  survey,  with  mean  and  standard  deviation 
(vertical  bars)  values  for  KB  and  ELS  sites;  (a)  Ocdom.  (b)  onap.  (c)  04,,  (d)  a,  by  summation  of 
Ox,  and  (e)  contributions  of  a^.  Note  scaling  difference  for  a^  and  a  in  (a)-(d). 
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Tables 

Values  of  a,  and  relationships  with  optically  active  constituents,  as  ratios,  for  2006  survey. 


acDOM(440) 

flAMp(440) 

aMp(440) 

TSS 

aMp(440) 

FSS 

ap(440) 

TSS 

□.^(440) 

a*  (676) 

0^(440) 

Chi 

0^(676) 

Chi 

(in-') 

(m-) 

(m2g-») 

(m2g->) 

{m2g-i) 

(m-') 

(m-«) 

(m2  mg-’) 

(m2  mg  ') 

KWO 

1.607 

0247 

0.1279 

0.2472 

0.2065 

0.152 

0.073 

0.0395 

0.0190 

KWl 

0.674 

0.113 

0.0791 

0.2833 

0.1325 

0.077 

0.026 

0.0470 

0.0161 

KB1 

0.180 

0.020 

0.0673 

0.1010 

0.1558 

0.027 

0.010 

0.0553 

0.0218 

KB2 

0267 

0.022 

0.0928 

0.1113 

0.1843 

0.022 

0.007 

0.0407 

0.0128 

KB3 

0.252 

0.027 

0.0942 

0.1232 

0.1602 

0.019 

0.004 

0.0445 

0.0105 

KB4 

0.186 

0.016 

0.0249 

0.0525 

0.0629 

0.024 

0.008 

0.0535 

0.0188 

KBS 

0218 

0.009 

0.0302 

0.0518 

0.0959 

0.020 

0.008 

0.0411 

0.0167 

KB6 

0.173 

0.011 

0.0308 

0.0572 

0.1064 

0.028 

0.013 

0.0653 

0.0299 

ELSl 

0.183 

0.010 

0.018 

0.007 

0.0263 

0.0103 

ELS2 

0.105 

0.0065 

0.029 

0.011 

0.0504 

0.0199 

ELS3 

0.155 

0.0104 

0.027 

0.013 

0.0331 

0.0161 

Mean 

0.186 

0.016 

0.0732 

0.106 

0.1441 

0.027 

0.010 

0.0445 

0.0167 

Std.  dev 

0.439 

0.073 

0.037 

0.089 

0.0479 

0.041 

0.020 

0.0109 

0.0056 

400  500  600  700 

Wavelength  (nm) 


Fig.  4.  Spectral  contributions  (%)  of  a*  too  for  selected  sites:  (a)  KWO.  (b)  KWl.  (c)  KB4. 
and  (d)  ELS2. 


squares  regression  (Fig.  5a;  p<0.01 ).  Inclusion  of  the  river  and  inflow 
zone  sites  (see  inset.  Fig.  5a)  results  in  a  stronger  relationship,  but  this 
is  an  artifact  of  the  disproportionate  effect  of  their  much  higher 
values.  Higher  Unap  per  unit  PAVm  is  indicated  for  these  KW  sites.  The 
relationship  between  Unap  ^nd  FeXRI  for  the  KB  and  ELS  sites  was 
somewhat  stronger  (Fig.  5b:  r^  =  0.82.  p<0.01). 

The  relationship  between  the  Mie  theory  estimates  of  am(440)  for 
n'  =  0.008  and  PAV^  is  presented  (Fig.  5c)  with  a  linear  least  squares 
regression  fit,  along  with  fits  for  am(440)  calculations  for  the  other  n' 
values  considered  (0.003.  0.006,  0.010).  The  estimates  of  Om  are 
shown  to  be  approximately  linearly  related  to  PAV^  for  all  samples 
(r^  =  0.97,  p<0.01)  and  when  KWO  and  KWl  are  not  considered 
(inset  of  Fig.  5c;  r^  =  0.62,  p  =  0.02).  Moreover,  increases  in  a  non¬ 
linear  manner  as  n'  increases  (Fig.  5c).  The  relationship  between  Onap 
(440)  and  the  0^(440)  estimates  is  depicted  for  the  scenario  of  n'  = 
0.008  (Fig.  5d).  According  to  this  scenario  am(440)  would  represent 
on  average  84%  of  aNAp(440).  About  40%  of  the  variability  in  aNAp(440) 
was  explained  by  the  independent  estimates  of  0^(440)  for  n'  =  0.008 
(p  =  0.10).  For  the  n'  =  0.006  case  (not  shown),  0^(440)  estimates 
represented  on  average  about  79%  of  aNAp(440).  Positive  y-intercept 
values  were  obtained  for  all  scenarios  that  converged  at  a  value  of 
about  0.006  m“  \  though  these  were  not  significantly  different  from 
zero. 

Discussion 

Spectral  characteristics 

The  slopes  of  the  CDOM  and  NAP  absorption  spectra,  Scdom  ^nd 
Snap,  generally  considered  to  be  proxies  of  the  composition  of 
contributing  constituents  (Babin  et  al.,  2003b).  Scdom  has  been  used  to 
correct  remote  sensing  bio-optical  data  (Schwarz  et  al..  2002). 
Changes  in  Scdom  h^ve  been  attributed  to  differences  in  CDOM 
composition,  sources  and  sinks  (Helms  et  al.,  2008),  and  changes  in 
molecular  weight  of  fulvic  acids  (Carder  et  al..  1989;  Blough  and 
Green,  1995).  Both  Increases  (Vodacek  et  al.,  1997;  Del  Vecchio  and 
Blough.  2004)  and  decreases  (Morris  and  Hargreaves.  1997)  in  Scdom 
from  photodegradation  have  been  reported. 

The  variety  of  protocols  adopted  to  estimate  Scdom  from  measured 
spectra  has  made  comparison  of  results  from  different  studies 
difficult,  as  the  effects  are  large  relative  to  the  real  differences  that 
may  be  attributable  to  differences  in  CDOM  composition  (Twardowski 
et  al..  2004;  Helms  et  al..  2008;  Loiselle  et  al..  2009).  The  two  primary 
sources  of  variability  associated  with  these  protocols  are  the 
wavelength  interval  considered  and  the  calculation  method  adopted 
(Twardowski  et  al.,  2004).  Narrow  wavelength  intervals  (e.g.,  100- 
150  nm  widths)  have  advantages  (Helms  et  al..  2008)  and  have 
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generally  been  adopted  in  recent  studies  (Babin  et  al..  2003b;  Binding 
et  al..  2008;  Perkins  et  al..  2009. 2010).  Some  modest  deviations  in  the 
wavelength  specifications  persist,  though  their  effects  are  probably 
minor;  e.g..  350-500  nm  by  Babin  et  al.  (2003a).  350-440  nm  by 
Binding  et  al.  (2008),  and  400-500  nm  by  Perkins  et  al.  (2009.  2010; 
and  adopted  here).  A  unification  of  calculation  protocols  appears  to 
have  evolved  recently  that  supports  the  adopted  non-linear  regres¬ 
sion  fitting  approach  (Babin  et  al.,  2003b;  Twardowski  et  al..  2004; 
Perkins  et  al..  2009.  2010). 

Most  of  the  Scdom  observations  for  the  bay  and  Lake  Superior 
(Tables  2  and  3)  were  at  or  slightly  below  the  lower  bounds  of  the 
populations  reported  elsewhere.  A  single  Ocdom  spectrum 
(ScDOM-O.Oll)  presented  by  Li  et  al.  (2004;  fig.  8.  p.  454)  from  a 
site  (Fig.  1)  in  west  central  Lake  Superior  (1999  collection)  is 
consistent  with  our  observations  (Tables  2  and  3).  The  mean  of  our 
bay  and  main  lake  observations  for  the  2006  survey  (0.0126  nm^ ') 
and  the  average  for  the  2007  survey  (0.01 1  nm“  ^)  were  lower  than 
those  reported  for  other  systems  to  date;  e.g.,  0.0176  nm“^  for 
European  coastal  waters  (Babin  et  al.,  2003b),  0.0161  nm"  ’  for  Lake 
Erie  (Binding  et  al..  2008),  0.0168  nm“  ’  for  Onondaga  Lake,  NY  (90 
observations.  Perkins  et  al..  2010),  0.0141-0.0186  nm~  ’  for  the  11 
Finger  Lakes  of  New  York  (Perkins  et  al.,  2010).  The  variability  among 
sites  encountered  for  this  attribute  in  this  study,  cv  =  14%.  was  similar 
to  that  reported  by  Babin  et  al.  (2003a)  (cv=11.4%)  for  European 
coastal  waters.  The  generally  low  values  of  Scdom  obtained  suggest  an 
atypical  composition  for  CDOM  in  this  system.  An  unusually  large 
contribution  of  high  molecular  weight  humic  acids,  relative  to  the 
lower  molecular  weight  fulvic  acids,  is  one  consistent  (Carder  et  al., 
1989)  potential  explanation.  Some  of  the  greatest  deviations  from  the 
mean  were  observed  for  near-shore  locations  (e.g.,  KBl,  Table  3;  SU22 
and  SU23  Table  4),  suggesting  local  differences  in  CDOM  composition. 
Binding  et  al.  (2008)  reported  even  greater  spatial  variations  in  Scdom 
in  Lake  Erie  based  on  90  observations. 

The  Scdom  value  for  the  humic-rich  Sturgeon  River  was  higher  than 
the  ELS  sites  and  most  of  the  KB  sites  (Table  3;  but  similar  to  five  other 
tributaries  to  Lake  Superior  sampled  in  April  2009,  unpublished  data. 
Upstate  Freshwater  Institute),  counter  to  expectations  based  on 
photoirradiation  effects  on  lake  CDOM  and  system-specific  slopes 
reported  by  Ma  and  Green  (2004).  albeit  for  a  different  wavelength 
range  (230-400  nm).  Their  river  value  (0.0144  nm  ^  ^ )  was  similar  to 
our  observation  (0.0169  nm~^).  However,  their  slope  value  for  lake 
water,  0.0259  nm~\  was  much  higher  than  reported  here  for  Scdom 
and  other  systems  where  longer  (visible)  wavelength  ranges  were 
considered,  suggesting  an  effect  of  the  different  wavelength  intervals 
considered. 

Snap  has  received  far  less  attention;  for  example,  values  have  been 
reported  in  only  a  limited  number  of  studies  of  lacustrine  systems 
(Binding  et  al.,  2008;  Perkins  et  al..  2009,  2010).  The  Snap  values 
reported  here  (mean  of  0.0129  nm~^)  are  similar  to  those  reported 
elsewhere,  where  the  same  wavelength  and  calculation  protocols 
were  used.  Babin  et  al,  (2003a)  reported  a  mean  of  0.0123  nm  ^ 
(cv=  1 1%).  Perkins  et  al.  (2009)  reported  an  overall  mean  of  0.0121 
for  the  Finger  Lakes  of  New  York,  with  means  for  the  individual  lakes 
in  the  range  of  0.01 13-0.0126  nm“  ^  The  mean  for  Onondaga  Lake 
was  0.0126  nm“  ^  (Perkins  et  al.  2010).  The  mean  for  Lake  Erie  was 
reported  to  be  somewhat  lower  (0.011  nm"\  Binding  et  al.,  2008). 
However,  the  wavelength  interval  and  calculation  protocols  were  not 
specified  for  the  Lake  Erie  observations.  The  Onap  spectra  reported 
here  had  shapes  (Fig.  2b)  that  were  also  generally  similar  to  those 
described  for  natural  minerogenic  particle  assemblages  that  have 
been  attributed  to  the  effects  of  iron  content  (Babin  and  Stramski, 
2004). 


Fig.  5.  Evaluations  of  relationships  related  to  minerogenic  particles  and  Qnap.  with  linear 
least-squares  regression  statistics:  (a)  Qnap  versus  PAVm.  (b)  Qnap  versus  average  FeXRI. 
(c)  Om  versus  PAVm,  and  (d)  On^p  versus  a^.  Insets  include  KWO  and  KWl,  except  for  (c). 
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Binding  et  al.  (2008)  presented  evidence  that  a  fraction  of  NAP  was 
bleached  by  the  sodium  hypochlorite  treatment  that  obscured  the 
440  nm  peak  in  Lake  Erie  samples  at  acDOM(440)  levels  as  low  as 
0.10  to  0.18  m~  \  The  effect  was  greater  as  Ocdom  increased  (Binding 
et  al.,  2008).  There  were  no  indications  that  this  was  a  problem  for 
resolution  of  a^,  spectra  for  Lake  Superior  and  Keweenaw  Bay,  as  the 
peaks  at  \  =  440  nm  were  well  resolved  in  all  samples  (e.g..  Fig.  2c). 
Moreover,  the  threshold  of  concern  for  acDOM(440)  for  the  effect  to  be 
manifested  (Binding  et  al.,  2008)  was  exceeded  in  most  of  the  2006 
Lake  Superior  samples  (Table  5).  The  values  of  the  ratio  a,j,(440) 
(676)  observed  for  the  KB  and  ELS  sites  were  well  centered  in  the 
portion  of  the  population  of  observations  reported  for  European 
coastal  waters  with  similar  levels  of  Chi,  with  the  exception  of  the 
highest  ratio  obtained  at  KBS  (Table  5).  The  high  ratio  for  this  sample 
is  probably  spurious  reflecting  the  diminished  accuracy  of  the  04,(676) 
peak  from  instrumentation  noise.  Variations  in  spectral  features, 
including  04,(440) 04,(676),  are  most  often  attributed  to  differences 
in  the  packaging  effects  and  accessory  pigment  content  of  phyto¬ 
plankton,  features  that  depend  on  community  composition  (Bricaud 
et  al..  1995;  Babin  et  al.,  2003b). 

patterns  and  relationships  with  optically  active  constituents 

The  general  spatial  patterns  of  Ocdom.  Qnap.  and  PAV^  in  the  2006 
survey  are  consistent  with  terrigenous  sources  of  these  constituents, 
though  the  details  are  doubtless  subject  to  substantial  variations.  For 
example,  Budd  (2004)  resolved  extensive  spring-time  sediment 
plumes  in  the  Keweenaw  region  of  Lake  Superior  through  satellite 
remote  sensing.  Cons  and  Auer  (2004)  reported  similar  variability  in 
flcDOM  (440)  (single  wavelength  measurements  only  within  the  bay) 
in  a  2003  survey,  though  their  values  were  shifted  somewhat  lower. 
The  results  of  the  2007  survey,  as  well  as  those  for  the  ELS  sites  in 
2006,  indicate  noteworthy  variations  in  Qcdom  also  occur  in  the  open 
waters  of  the  lake.  Binding  et  al.  (2008)  reported  even  greater  spatial 
variability  in  Lake  Erie,  with  a  range  of 0.08-0.75  m”  ^  for  acDOM(440), 
compared  to  0.105-0.267  m~  ^  for  Keweenaw  Bay  and  Lake  Superior 
sites  reported  here.  The  mean  for  the  bay  and  Lake  Superior  sites 
observed  in  this  study  was  25%  less  than  reported  for  Lake  Erie 
(Binding  et  al.,  2008).  The  partitioning  of  ax(440)  based  on  the  single 
set  of  spectra  presented  by  Li  et  al.  (2004;  fig.  8.  p.454)  tracked  closely 
that  reported  for  the  ELS  sites  (Table  5). 

Comparisons  of  specific  absorption  coefficient  values  obtained  for 
Ox  here  with  the  limited  observations  reported  for  other  systems 
establish  substantial  differences,  indicating  the  need  for  system- 
specific  relationships  to  support  mechanistic  modeling  goals.  The 
average  specific  absorption  coefficient  for  CDOM  based  on  DOC 
(a*cDOM.  Doc( 440)  =  0-107  m^g“^)  from  the  2007  cruise  was  similar 
to  the  mean  reported  for  Lake  Erie  (0.102  m^g~\  Binding  et  al., 
2008).  However,  the  range  for  Lake  Superior  (0.052-0.167  m^g'^) 
was  narrower  (0.035-0.306  m^  g~  \  Binding  et  al.,  2008).  The  utility 
of  DOC  as  an  optically  active  constituent  to  which  Ucdom  can  be  linked 
is  limited.  In  particular,  the  fraction  of  DOC  that  absorbs  light  (CDOM) 
is  subject  to  substantial  differences  between  systems  and  temporal 
variations  within  individual  systems  (Morris  and  Hargreaves,  1997). 
flcDOM  and  DOC  were  not  significantly  correlated  for  the  paired 
observations  of  the  2007  survey,  a  common  finding  (e.g.,  Perkins  et  al., 
2010).  The  a*cDOM.Doc(440)  was  nearly  2-fold  greater  in  Onondaga 
Lake,  NY  (Perkins  et  al.,  2010).  DOC  was  not  included  in  the 
comprehensive  survey  of  Ox  in  coastal  waters  around  Europe  (Babin 
et  al.,  2003b),  likely  in  recognition  of  the  disconnection  between  this 
aggregate  measure  and  Qcdom- 

Specific  absorption  coefficient  information  for  NAP  for  lacustrine 
systems  is  rare.  The  a*NAP.T5s(440)  and  a*NAP.Fss(440)  values  reported 
here  are  consistent  with  values  reported  elsewhere  for  cases  of 
dominance  of  Qnap  by  minerogenic  particles  in  marine  systems  (Babin 
et  al.,  2003b;  Babin  and  Stramski,  2004).  The  mean  a*NAP.TSs(440) 


value  for  the  KB  sites  (0.057  m^  g^  ^)  was  similar  to  that  reported  for 
Lake  Erie  (0.0403  m^  g~  \  Binding  et  al.,  2008)  and  European  coastal 
waters  (0.041  m^  g“  \  Babin  et  al.,  2003b),  but  less  than  for  Onondaga 
Lake  (0.090  m^g~\  Perkins  et  al..  2010).  The  mean  a*NAP.FSs(440) 
value  for  the  bay  (0.083  m^  g~^)  was  -35%  lower  than  reported  for 
Lake  Erie  (Binding  et  al.,  2008),  but  was  2-fold  lower  than  for 
Onondaga  Lake  (Perkins  et  al.,  2010).  The  findings  of  Babin  and 
Stramski  (2004)  and  the  FeXRl  dependence  reported  here  (Fig.  5b) 
suggest  at  least  a  portion  of  the  differences  in  the  specific  absorption 
coefficient  for  NAP  between  studies  may  be  due  to  differences  in  Fe 
content  of  suspended  solids.  The  mean  a’'p.Tss(440)  for  the  bay 
(0.128  m^g~^)  was  -45%  lower  than  reported  for  Onondaga  Lake 
(Perkins  et  al..  2010). 

The  average  Chl-specific  a<p(440)  for  the  KB  and  ELS  sites,  a*4,.chi 
(440)  =  0.046  m^  mg~’,  was  generally  consistent  with  most  of  the 
reported  values  in  both  the  marine  (Bricaud  et  al..  1995;  Babin  et  al., 
2003b)  and  lacustrine  literature  (Bukata  et  al.,  1991;  Effier,  1996; 
Perkins  et  al..  2010).  However,  it  was  distinctly  lower  than  recently 
reported  for  Lake  Erie  (0.086  m^  mg Binding  et  al.,  2008).  The 
average  value  for  a*4,,chi(676),  0.017m^mg~\  was  similar  to  that 
reported  for  Onondaga  Lake  (Perkins  et  al..  2010). 

Minerogenic  particles  and  amp 

Empirical  analyses 

The  strong  relationship  between  Unap  and  PAV^  for  the  KB  and  ELS 
sites  (Fig.  5a)  supports  the  position  that  terrigenous  (Table  2) 
minerogenic  particles,  mostly  clay  minerals,  were  important  in 
regulating  the  observed  differences  in  that  component  of  a.  PAVm 
has  both  system-specific  and  conceptual  advantages  over  the 
gravimetric  measure  of  minerogenic  particles  (FSS).  PAV^  is  a  more 
sensitive  and  precise  (Peng  and  Effier,  2007;  Peng  et  al.,  2009)  metric 
than  FSS  in  dilute  systems  such  as  Lake  Superior.  Moreover,  particle 
projected  area  (e.g.,  PAV^)  is  more  fundamentally  connected  to  the 
light  attenuating  properties  of  particles  (Owen.  1974;  Iturriaga  and 
Siegel,  1989;  Peng  and  Effier,  2007)  than  gravimetric  representations. 

The  positive  y-intercept  (0.003  m  ^)  of  the  QNAP-PAVm  relation¬ 
ship  is  consistent  with  contributions  from  another  component(s)  to 
Qnap.  that  is  likely  organic  detritus  (Babin  et  al.,  2003b;  Binding  et  al., 
2008).  fliMAP  is  expected  to  be  well  represented  by  a  two  component 
model 

^^NAP  =  (b) 

where  Od  =  absorption  associated  with  particulate  detritus. 

The  even  stronger  relationship  between  Qnap  ^nd  FeXRl  (Fig.  5b) 
further  supports  the  position  that  minerogenic  particles  were 
important  in  regulating  the  observed  spatial  differences  in  Unap 
among  the  KB  and  ELS  sites.  Using  pure  minerogenic  specimens  and 
an  array  of  land-based  natural  deposits,  Babin  and  Stramski  (2004) 
reported  Fe  content  was  a  primary  driver  of  0^.  To  the  authors’ 
knowledge,  the  results  presented  here  for  Lake  Superior  are  the  first 
support  (Fig.  5b)  for  a  potentially  central  role  for  Fe  content  of 
minerogenic  particles  in  regulation  of  Qnap  based  on  analysis  of  water 
column  samples.  However,  the  possibility  that  the  observed  relation¬ 
ship  is  instead  an  artifact  of  a  covarying  factor  cannot  be  discounted. 
Absorption  by  Fe  has  been  associated  with  a  combination  of  the  effects 
of  Fe  adsorbed  to  particle  surfaces  and  embedded  as  an  impurity  in 
crystal  lattices,  with  spectral  absorption  shapes  similar  to  a  pure  iron 
hydroxide  suspension  (Babin  and  Stramski.  2004).  The  potential  for  a 
non-linear  relationship  was  indicated  for  our  results  by  the  negative 
y-intercept  for  the  linear  regression  expression  and  the  much  higher 
Onap  values  for  modest  increases  in  FeXRl  for  the  KWO  and  KWl 
samples  (see  inset.  Fig.  5b).  However,  higher  detritus  concentrations 
at  these  two  sites  may  also  have  contributed  to  the  latter  feature. 
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Uncertainty  prevails  with  regards  to  the  interplay  between  CDOM 
and  absorption  by  particles  in  general  (Binding  et  al.,  2008),  and 
minerogenic  particles  specifically  (Gallegos  and  Neale,  2002).  The 
extent  to  which  this  can  modify  the  inherent  absorption  character¬ 
istics  of  mineral  particles  (e.g.,  increase  the  effective  n')  remains 
unresolved. 

Mie  theory  estimates  of 

The  strong  relationships  between  Mie  theory  estimates  of  and 
measurements  of  PAVm,  illustrated  for  the  case  n'  =  0.008  (Fig.  5c), 
supports  the  position  that  PAV^  can  also  serve  as  a  surrogate  metric 
for  Om,  as  well  as  for  (Peng  and  Effler,  2007).  The  variability  in  the 
relationship  between  0^(440)  and  PAVm  is  attributable  to  differences 
in  the  particles  size  dependencies  of  these  features,  as  well  as  modest 
variations  in  the  contributions  of  the  various  chemical  classes  to  the 
minerogenic  particle  populations  of  these  samples  (Table  2).  Based  on 
the  limited  literature  (Patterson  et  al..  1977;  Wozniak  and  Stramski, 
2004),  the  scenarios  of  n'  =  0.006  and  0.008  appear  to  be  the  most 
reasonable  for  the  clay  particles.  For  example,  an  n'  value  of 
approximately  0.007  was  specified  for  A  =  440nm  in  a  spectrum 
(350-750  nm)  adopted  by  Wozniak  and  Stramski  (2004)  in  modeling 
the  optical  properties  of  mineral  particles  in  seawater.  Specification  of 
much  lower  values  for  quartz  and  Ca-rich  (e.g.,  calcite)  particles,  as 
adopted  here  (Table  2),  is  consistent  with  Um  measurements  on  pure 
samples  (Babin  and  Stramski,  2004). 

A  rigorous  test  of  optical  closure  of  the  forward  estimates  of  Om 
(440),  based  on  integration  of  individual  particle  information  from 
SAX  into  Mie  theory  calculations,  cannot  be  supported  by  this  data  set 
because  Oa  was  not  independently  quantified  (Eq.  (2)).  However,  the 
relationship  between  aNAp(440)  and  these  forward  method  estimates 
of  am(440),  presented  here  (Fig.  5d)  for  the  first  time,  depicts  features 
of  consistency  that  are  encouraging  for  future  SAX-based  closure 
efforts  and  partitioning  of  Onap  according  to  Om  and  Od.  The  level  of 
consistency  is  impressive  in  the  context  of  the  multiple  sources  of 
uncertainty  that  could  influence  the  results.  First,  the  estimated 
magnitudes  of  am(440),  that  incorporated  reasonable  value(s)  of  n', 
were  generally  consistent  (e.g.,  appropriate  order  of  magnitude)  with 
the  measurements  of  aNAp(440)  (Fig.  5d).  Secondly,  the  best  fit  linear 
least  squares  regression  relationship  (r^  =  0.40)  approached  signifi¬ 
cance  (p  =  0.1).  and  was  generally  consistent  with  indications  from 
the  empirical  PAVm-based  model  (Fig.  5a).  Third,  the  indicated 
positive  y-intercept  is  consistent  with  a  modest  Od  component  that, 
while  poorly  quantified,  is  generally  acknowledged  to  be  present, 
even  in  case  1  waters  (Iturriaga  and  Siegel,  1989;  Babin  et  al.,  2003b; 
Wozniak  and  Stramski.  2004). 

Summary  and  implications 

This  study  has  provided  the  most  complete  characterization  of  a 
and  Qx  in  Lake  Superior,  confirming  that  it  is  an  optically  complex  case 
2  system.  The  findings  have  extended  the  limited  literature  for  Ox  and 
relationships  with  optically  active  constituents  for  lacustrine  systems 
and  particularly  for  the  Laurentian  Great  Lakes.  Differences  among 
systems  support  the  continued  need  for  system-specific  characteriza¬ 
tions  as  presented  here.  The  spectral  features  of  Unap  (Snap)  ^nd  (04, 
(440)^04,(676))  were  similar  to  those  reported  for  coastal  marine 
and  lacustrine  case  2  systems.  However,  the  Scdom  values  for 
Keewenaw  Bay  and  open  waters  of  Lake  Superior  were  low  relative 
to  observations  reported  for  other  case  2  waters  (and  similar  to  Snap). 
suggesting  an  atypical  composition  for  the  lake's  CDOM.  Ocdom  was 
the  dominant  absorbing  component  through  the  blue  and  green 
wavelengths,  representing  >75%  of  a(440)  at  all  monitored  locations 
in  2006.  Spatial  patterns  of  Ocdom  and  Qnap  indicated  the  effects  of 
terrigenous  sources.  A  general  gradient  in  Qcdom  and  Unap  extended 
from  the  Sturgeon  River,  through  the  bay,  and  into  eastern  Lake 
Superior  in  the  2006  survey.  In  contrast,  04,  was  similar  in  Lake 


Superior  and  the  bay  during  the  survey  of  2006.  Differences  in  Fe 
content  (e.g.,  FeXRI  ^  PAVm)  probably  contributed  to  variations  in  the 
relationship  between  Onap  and  suspended  solids,  as  represented  by 
a*NAP.Tss(440)  and  a*NAP.Fss(440),  reported  here  and  elsewhere. 

The  absorption  characteristics  of  Lake  Superior  (Figs.  2-4)  present 
difficult  challenges  for  effective  remote  sensing  of  Chi  patterns,  that 
are  shared  by  many  other  lacustrine  oligotrophic  case  2  waters.  The 
blue  absorption  peak  is  obscured  by  the  other  absorbing  components, 
particularly  CDOM  (e.g..  Fig.  4c  and  d).  The  red  (\  =  676)  peak 
alternative  to  estimate  Chi  may  be  compromised  because  this  signal  is 
overwhelmed  by  (e.g..  Fig.  4c  and  d).  These  conditions  have 
doubtless  contributed  to  the  reported  problems  in  remote  sensing  of 
Chi  in  the  lake  (Budd  and  Warrington,  2004;  Li  et  al..  2004).  The 
similarity  of  Scdom  and  Snap  make  it  difficult  to  resolve  the 
contributions  of  these  two  components  through  remote  sensing. 

The  characterizations  of  Ox.  and  the  relationships  with  optically 
active  constituents,  presented  here,  together  with  the  recently 
reported  backscattering  information  (Peng  et  al.,  2009),  will  serve 
to  support  more  robust  remote  sensing  algorithm  development  and 
advance  mechanistic  approaches  for  Lake  Superior.  These  character¬ 
izations  meet  specific  recommendations  made  to  support  remote 
sensing  initiatives  for  the  lake  (Gons  and  Auer,  2004;  Gons  et  al.. 
2008).  Moreover,  this  information  can  support  parameterization  of  a 
mechanistic  optics  model  for  the  underwater  light  field  to  simulate 
AOPs  such  as  Secchi  depth  and  the  diffuse  attenuation  coefficient  for 
downwelling  irradiance  (e.g.,  Effler  et  al..  2008). 

This  study  has  advanced  the  understanding  of  the  role  of 
minerogenic  particles  in  contributing  to  Qnap.  based  on  analysis  of 
paired  SAX  results.  A  widely  held  notion  has  been  that  Unap  is 
dominated  by  in  case  1  waters  as  well  as  in  many  case  2  systems 
(Morel  and  Ahn,  1991;  Babin  et  al.,  2003b).  The  strong  relationships 
between  Qnap  ^nd  PAV^  (Fig.  5a),  and  Qnap  and  FeXRI  (Fig.  5b), 
reported  here,  provide  compelling  empirical  evidence  for  the  central 
role  of  terrigenous  minerogenic  particles  in  regulating  the  observed 
differences  in  Qnap  amongst  these  sites.  In  particular,  this  FeXRI 
relationship  represents  its  first  demonstration  based  on  analysis  of 
water  column  samples.  Heretofore,  such  a  relationship  had  only  been 
anticipated,  based  on  analysis  of  land-based  sediment  samples  (Babin 
and  Stramski.  2004).  A  central  role  of  minerogenic  particles  was 
further  supported  by  the  consistency  of  predicted  Om  values  with  Unap 
observations,  obtained  from  the  first  set  of  forward  method  Mie 
theory  calculations  based  on  SAX  data  and  reasonable  literature 
values  for  n'. 
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